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ABSTRACT 

Context. Although most young massive stars appear to be part of multiple systems, it is poorly understood how this multiplicity 
influences the formation of massive stars. The high-mass star-forming region W3 IRS5 is a prime example of a young massive cluster 
where the cluster center is resolved into multiple subsources at cm and infrared wavelengths, a potential proto-Trapezium system. 
Aims. We investigate the protostellar content in the 1.4 mm continuum down to subarcsecond scales and study the compact outflow 
components, also tracing the outflows back to their driving sources via the shocktracing SiO and SO2 emission. 
Methods. The region W3IRS5 was mapped with the PdBI at 1.4 mm and 3.4 mm in the AB configurations, tuning the receivers to 
observe the molecular transitions SO2 (222_2o-22i,2i), SO2 (83,5-92,8), SiO (2-1), and SiO (5-4). 

Results. In the continuum we detect five sources, one of them for the first time, while counterparts were detected in the NIR, MIR 
or at radio wavelengths for the remaining four sources. Three of the detected sources are within the inner 2100 AU, where the 
protostellar number density exceeds 10* protostars pc"^ assuming spherical symmetry. Lower limits for the circumstellar masses of 
the detected sources were calculated, ranging from ~0.3 to ~40 Mq although they were strongly affected by the spatial filtering of the 
interferometer, losing up to ~90% of the single-dish flux. However, the projected separations of the sources ranging between ~ 750 
and ~ 4700 AU indicate a multiple. Trapezium-like system. We disentangled the compact outflow component of W3 IRS5, detecting 
five molecular outflows in SiO, two of them nearly in the line of sight direction, which allowed us to see the collapsing protostars 
in the NIR through the cavities carved by the outflows. The SO2 velocity structure indicates a rotating, bound system, and we find 
tentative signatures of converging flows as predicted by the gravoturbulent star formation and converging flow theories. 
Conclusions. The obtained data strongly indicate that the clustered environment has a major influence on the formation of high-mass 
stars; however, our data do not clearly allow us to distinguish whether the ongoing star-forming process follows a monolithic collapse 
or a competitive accretion mechanism. 

Key words, stars: formation - instrumentation: high angular resolution - instrumentation: interferometers - ISM: individual objects: 
W3 IRS5 

1. Introduction quence. The radius where these Hii regions become hydrody- 

^, ^ ■ r, ■ , namically supported -thus stopping the accretion flow- scales in- 

The formation of high-mass stars appears to be intimately hnked ^^^^^j ^.^^ attracting mass, thus massive trapezia systems 
^ with the formation of multiple stellar systems. Optically vis- allow more m ass to be accreted to the central accreting objects 

■ ible OB -type stars show a much highe r degree of multiplic - ( |Ketoll200"2allbl 120^1) 
" ' ' ity than their lower-mass counterparts dPreibisch et al.l Il999h . j 1 ■ j • 1 

. / , ■ • , , — 7^, 7^ Second, binary systems are expected to form in the center 

Now observations are showing that such high-mass multi- , \ rr • t-i. i.- • n 

, , . . , ° „ „ „ „ . , , , of these dense proto-Trapezia. These binaries will continue ac- 

ple systems have separations ranging between 10,000AU and . r 1 j . . i. i.- j . 1 

V.TT J 11 J i-i 1 creting gas, funneled to the center by the combined potential 

1 AU, and appear to be bound groups within larger clusters ^ , , r, , , r , ■ 

rrj . 0^ o 1- , r , of the protocluster. Because the angular momentum of the in- 

OMerrmlliod & Garciall2001D . In contrast to hierarchical svstems , ,. ii ^c.i- 

, : — . , . , . falling gas has no correlation with that of the binary s, the masses 

where the successive separation between Its members increases r,. • j- j 1 . • ui l- r i.- 

, , „ , J • 11 of the individual stars increase while the separation of the bi- 

by large factors, many of these systems are dynamically un- , . . , 11 1 

, , , . 1 • 1 1 r , 11 J nary decreases, giving as a natural outcome a high-mass close 

stable, nonhierarchical clusters of three or more stars, called , . 1 1 j n j 

. „ , „ . 1 ^ ■ ^ ■ r,,, , , binary system that could eventually undergo mergers, creating 

trapezia after the Trapezium cluster in Orion. The role such sys- . ^ ,„ 11 o n f ^/^rvcT 

, . ■ r ■ ■ J j i even more massive systems (Bonnell & Bate 2005). 

tems play in massive star formation is not yet understood, al- 

though several suggestions have been made. . ^he previous examples show how the formation of mas- 

tr ^. . 1 » u u- u sive Stars in multiple systems with separations on the order of 

First, the mass density of these multiple systems may be high , , rr ,, .• • ■. , 

■ . 11 t u * 11 1 ~1000AU may be affected by their companions, either by mod- 
enough to gravitationally trap hypercompact or even small ul- . ^ . , , , ^ . ^ , . 
tracompact Hn regions. S & has shown that ^^ymg the morphology of accreting envelopes allowing more 

.. . »u u »u ui^u ■ u • mass to be accreted, or by forming more massive stars by merg- 

continued accretion through the HCHii region may be possi- . , . , , ° .7 & 

ui f» »u ? » u u J »u • ing high-mass close binaries, 

ble, even after the supporting stars have reached the main se- ° ° 

Some multiple high-mass protostellar systems are known. 

* FITS images are available in electronic form at the CDS via For ex ample, W 33A was imaged by van der Tak & Mentei 

anonymousftptocdsarc.u-strasbg.fr (130.79.128.5) or via (l2005|) at 43 GHz with the Very Large Array findii ig three con- 

|http : //cdsweb . u-strasbg . fr/cgi-bin/qcat?3/A+A/, tinuum sources at separations of 3000 - 5000 AU; iHunter et al.l 
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( I2OO6I) found similar systems in NGC 6334 1 and NGC 6334 1(N) 
with a few millimetric sources inside a region of ~1 0,000 AU. 
More recently, Beutheretal. (2007) resolved the central 
7,800 AU of the hot molecular core G29.96-0.02 into four sub- 
millimeter peaks, a potential proto-Trapezium system. 

One prime example of a trapezia is the system W3 IRS5 in 
the W3-Main region of the Hii/molecular cloud complex W3, 
a very active star-forming region part of the larger star-forming 
complex W3/W4 located in the Perseus arm of the Galaxy, that 
also encompasses the W4 star-forming region and the OB asso- 
ciations IC 1805 and IC 1795. In W3 the star formation takes 
place in the embedded regions W3-Main, W3-North, and W3- 
OH 

lOev et akl (|2005|) found evidence supporting the scenario in 
which the star formation in W3-Main was triggered by the for- 
mation of IC 1795, which in turn was triggered by the Perseus 
chimney/superbubble. This scenario, however, has been brought 
into question by the recent Chandra X-Ray observations of 
iFeigelson & Townslevl (2008), because W3-Main does not show 
the elongated and patchy structure of a triggered star cluster 
W3-Main itself may harbor its own triggered star-formation sce- 
nario because the various morphological classes of Hn regions 
that can be found in it suggest that different stages of massive 
star formation (MSF) may be seq uentially triggered by t he pre s- 
sure of the expanding Hn regions (iTieftrunk et al.lfT997l fT998l) . 

The system W3 IRS5 is a known dou ble infrared (IR) source 
dHowell et al.1 119811; iNeugebauer et all IT9 82) with a total lu- 
minosity of 2 X 10^ (ICampbell et al.l "n995). located in the 
W3-Main region at a dis t ance o f 2kpc iMegeath et al. 2008). 
Recently, iMegeath et al.l (l2005l) presented HST observations 
with an angular resolution of 350 AU identifying seven near- 
IR sources including the two previously known. Three of these 
sources have both mid-IR c ounterparts a nd are 1.2 and 0.7 cm 
continuum sources jWilson et al. 2003; v an der Tak et al. | l2005l) . 
The sizes of the Hii regions are <240 AU, indicating that th ey are 
proba bly gravitationally bound and may be accreting (e.g.. lKetol 
I2003h . Based on its size, the probable ma sses of the stars, and 
its nonhierarchical distribution of sources, Megeath et al.l (l2005h 
proposed that this system is a proto-Trapezium system, which 
would emerge as a bound Trapezium similar to that in the Orion 
nebula. T his region is also the sourc e of at least two outflows 
dlmai et al. 2000; Wils on et all l2003t). and is in the center of 
an em bedded cluster of 80-240 low-mass stars dMegeath et al.l 
Il996t) which in turn is surrounded by a core of several hundred 
low-mass stars dFeigelson & Townslevll2008l) . 



2. Observations 

We have mapped W3IRS5 with the IRAM Plateau de Bure 
Interferometer (PdBlf] in the A (Jan/Feb 2006) and B (Mar 
2006) configurations at 1.4 and 3.4 mm imaging the contin- 
uum, achieving an angular resolution with a synthesized beam 
of 0. 39" X 0. 34" at 1 .4 mm and 1 .00" x 0. 89" at 3 .4 mm. At the 
given distance of 2 kpc that translates to a spacial resolution of 
-700 and -1900 AU, respectively. The ~ 0.36" resolution ob- 
tained was feasible because of the new very extended baselines, 
extending to over 700 m. 

The receivers were tuned to 87 and 217 GHz both in the 
lower sideband. With this spectral setup we observed the SO2, 
SiO and C'^O transitions described in Table[T]with a maximum 



spectral resolution of 0.5 km s '.We adopted a systemic veloc- 
ity Vlsr = -39kms-i ('Ridge & M oorell2001l: ITieftrunk et all 
1998). Although observed, C'**0 was not detected due to the 
short-spacings problem inherent to the interferometer filtering 
out the received flux. 



Table 1. Observed molecular transitions and rms of the respec- 
tive maps. 



Transition 


y 

(GHz) 


Spect. Resol. 
(MHz) 


(K) 


rms^ 
(mJy beam 


SO2 (83,5-92,8) 


86.64 


0.14 


55 


8.0 


SiO (2-1) 


86.85 


0.14 


6 


7.0 


Continuum 


87.02 






1.8 


SO, (222,2o-22,,2,) 


216.64 


0.36 


248 


18.0 


Continuum 


216.87 






1.2 


SiO (5-4) 


217.10 


0.36 


31 


14.0 


C"*0(2-l]^ 


219.56 
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' IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



" For a spectral resolution of 0.5 km s ' except for SiO (5-^), whose 
map has a spectral resolution of 1 km s"' . 
* Not detected, filtered out by the interferometer. 



The phase and amplitude calibrators were 0355+508 and 
0059+581 and the flux calibrator was 3C345, adopting the val- 
ues from the SMA flux monitoring of these quasars. The data 
were calibrated with the program CLIC and then imaged with 
the program MAPPING, both part of the GILDAS package. The 
spectra were processed with the program CLASS77, also from 
the GILDAS package. 

After inversion and cleaning, the continuum data have rms 
noise levels cr - 1.8 mJy beam"' and cr - 1.2mJybeam"' 
at 3.4mm and 1.4mm respectively. For the line data, with 
a spectral resolution of 0.5 km s ' the noise levels are 
cr - 7mJybeam"' for SiO (2-1), cr - 18mJybeam"' for 
SO2 (222,20-22i,2i) and cr = 8mJybeam-' for SO2 (83,5-92,8), 
while SiO (5^) was imaged with a spectral resolution of 
1 kms"' resulting in a rms cr - 14mJybeam"' (see Table[T]). 



3. Results 

3.1. Large-scale Overview 

The W3 complex is a very rich one when it comes to Hii re- 
gio ns. Several Hii region s have been detected at 6, 2 and 1.3 cm 
by ' Tieftrunk et al.l (II997 ). encompassing hypercompact, ultra- 
compact, compact and diffuse Hii regions. All of them probably 
originated from the same volume of neutral gas. 

With the ages they estimated, ITieftrunk et al.l d 19971) sug- 
gest a qualitative evolutionary sequence of the Hn regions, going 
from hypercompact to ultracompact, compact and then diffuse. 
Also the number, age and spatial distribution of the Hii regions 
suggest a scenario of triggered star formation by the expansion 
of the Hi: regions that formed first. 

In our 3.4 mm map we detected three of those Hii regions, the 
compact W3 B, the ultracompact W3 F and the hypercompact 
W3 IRS5, also known as W3 M, the subject of our study (see 
Fig. [TJ. The shell-like structure and overall morphology of W3 
B is recovered by our observations, as well as the cometary shape 
ofW3F 
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R.A. 


Dec. 


T " 


Sv" 


Mass (Mo) 




N(H2t 


Av 


MM 


(J2000) 


(J2000) 


(mJy beam"') 


(mJy) 


50 K 100 K 


200 K 


(10'" cm-2) 


(10' mag) 


1.4 mm with 0.39" x 0.34" beam 


1 ... 


02 25 40.77 


62 05 52.49 


34.8 


70 


6.7 3.2 


1.6 


8.0 


8.5 


2 ... 


02 25 40.68 


62 05 51.53 


20.6 


24 


2.3 1.1 


0.5 


4.7 


5 


3 ... 


02 25 40.66 


62 05 51.95 


10.4 


13 


1.2 0.6 


0.3 


2.4 


2.5 


4 ... 


02 25 40.75 


62 05 50.45 


5-71 


5.7^ 


0.61 0.3" 


0.1" 


1.3 


1.4 


3.4 mm with 1.00" x 0.89" beam 


1 ... 


02 25 40.77 


62 05 52.49 


23.1 


24 


83.4 40.1 


20.2 


30 


32 




02 25 40.68 


62 05 51.85 


11.2 


11.2'' 


39.0" 19.1" 


9.4" 


14 


15 


6 ... 


02 25 40.73 


62 05 49.86 


5.1 


5.1'' 


17.7" 8.7" 


6.3" 


6.6 


7 



" The fluxes have been corrected for the free-free contribution, except for source MM4 

* Calculated assuming T = lOOK. 

" Free-free contribution not subtracted. 

'' Unresolved, ly given instead of S,,. 

" Calculated with the peak intensity instead of the flux density, because the angular size of the source is smaller than the beam. 
^ Source MM5 is likely the merged pair MM2 and MM3. 



3.2. Millimetric Continuum 

Figure|2]shows the 1 .4 and 3.4 mm continuum maps, with the six 
sources detected labeled MMl-6 in descending order of inten- 
sity at 1.4mm for and at 3.4mm for MM5 and MM6. 
MMl is the only source resolved at both wavelengths, while 
MM4 is detected only at 1.4mm and MM6 only at 3.4mm. 
At the long wavelength we do not resolve MM2 nor MM3. 
However, based on its position, MM5 is likely the joint contri- 
bution from MM2 and MM3 and is not considered an individual 
source. 

The properties of the sources are summarized in Table |2] 
Columns 2 and 3 give their positions; the measured peak flux 
intensity and integrated flux density, both corrected for the free- 
free contribution, are given in columns 4 and 5 respectively. 

In the standard model of optically-thick (spectral index ~ 2) 
and optically-thin (spectral in dex ~ -0.1) free-free e mission, 
the spectral indexes derived bv lvan der Tak et a n (l2005h indicate 
that the 43 GHz emission is close to the "turnover point" be- 
tween optically-thick and optically-thin free-free emission, and 
is approximately an upper limit to the free-free contribution at 
millimetric wa velengths. In HCHii region r nodels with density 
gradients (e.g.. lKetoll2003l) or dumpiness (ISewilo et al.ll2004b 
the free-free emission may rise somewhat more, but unlikely by 
a significant amount. Therefore, here we assume that the 43 GHz 
emission is a good measure of the free-free contribution at milli- 
metric wavelengths and we use it to correct our measured fluxes 
for such contribution, as shown in Table [3] T his could be done 
for aU our sources except MM4, for which Ivan der Tak et al.l 
(I2005h did not detect a counterpart. 

Because the brightness temperature of the coiTesponding 
Planck function for the strongest source -MMl- is about 1 1 K, 
just ~IQ% of the typical hot core temperatures of ~100K, 
we can assume that the emission comes from optically thin 
dust and thus calculate th e masses an d colum n densities with 
the approach outlined bv [Hildebrandl (Il983h and adapted by 
[Beuther et al. (2002a, 20051) . We adopted a distance of 2kpc 
(Megeat h et al.il200 8) and used a grain emissivity index yS = 2, 
corresponding to a dust opacity per unit mass KiA,„„, ~ 0.05 and 
KiAmm ~ 0.3cm^g"' for a median grain size a - 0.1 //m and 
grain mass density p = 3 g cm"^. In columns 6-8 and 9 of Table 



Table 3. Contribution of the free-free emission to the measured 
millimetric fluxes. 



Peak Intensity (f-f) Flux Density (f-f) 

MM (mJybeam-') I' l I (mJy) 

' ^ measured ly | ^ ' measured Sy 

1.4 mm with 0.39" x 0.34" beam 



1 ... 


2 0.05 


2 


0.03 


2 ... 


2 0.09 


4 


0.14 


3 ... 


2 0.16 


3 


0.19 




3.4 mm with 1.00" x 0.8 


9" beam 




1 ... 


2 0.08 


2 


0.08 


5 ... 


4 0.26 


4* 


0.26'' 


6 ... 


0.7 0.12 


0.7'' 


0.12'' 



" Estimated fromlyan derTaketal] ( l2005h . see discussion in the text. 
Unresolved in our data, free-free emission correction for ly used. 



|2] are the calculated masses and H2 column densities respec- 
tively. We did not derive a temperature for this region, therefore 
the masses were calculated assuming the different temperatures 
T ~ 50K, T ~ lOOK and T ~ 200 K, while for the N(H2) calcu- 
lation we assume a hot core temperature of ~ lOOK. 

The obtained masses are strongly affected by the spatial 
filtering in herent to interferometer s. Comparing with single- 
dish data dvan der Tak et alj l2000l) we recover only ~ 9% 
of the flux at 1.4 mm (see Sec. 14.11 for further discussion). 
The visual extinctions in column 10 were calculated assuming 
Av = N (H2) /0.94 X 10^1 (Frerking et al...l982) . 

3.3. Silicon Monoxide emission 

Silicon Monoxide (SiO) is believed to trace strong shocks 
in dense molecular gas, particularly in molecular outflows 
(ISchilkeetal.lll997h . Here two SiO transitions were observed, 
SiO (5^) with the 1.4 mm receiver and SiO (2-1) with the 
3.4 mm receiver, seen in Figure [3] The five molecular outflows 
we detect are also marked, labeled SlO-a - e, four of them at 
both wavelengths and the larger, westernmost outflow, only at 
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3.4 mm. In Figure |4] we present the spectra of the five outflows, 
taken toward the center of the outflow for SlO-a, -b and -d, and 
toward the wings for SIO-c and -e. From the integrated spectra 
of the 1 .4 mm emission within a 4" x 4" region centered at the 
phase center (Fig.|4f), we determine that the contribution of the 
ambient gas emission ranges from -43 to -36 km s"' . This is also 
in agreement with the linewidth expected from the ambient gas 
emission for dense and small cores (Garay & Lizano 1999). 

The spatial coincidence of the blue- and redshifted emis- 
sion of SlO-a with MMl indicates that it is aligned very 
close to the line of sight (l.o.s.). This can be clearly seen in 
the SiO (5-4) transition at 1 .4 mm, where the most compact 
emission is mapped (Fig. [3j. At this wavelength the red- and 
blueshifted contours are overlapped, however in the SiO (2-1) 
transition at 3.4mm only a blueshifted emission peak is seen. 
There is overlapping redshifted emission, but it does not peak 
at the same position but rather ~0.5" to the southeast of the 
blueshifted peak. The position of the SiO (5^) peaks as well 
as the blueshifted SiO (2-1) peak match the position of MMl, 
and this is likely the driving source of this outflow. The spec- 
trum at the center position (Fig. |4^) exhibits a broad emission 
with a FWHM linewidth of 7 km s ^ showing also a strong peak 
at the systemic velocity and a redshifted narrow (FWHM of 1.3 
km s"' ) feature peaking at ~ -32 km s"' . 

For SlO-b on the other hand, we detect and resolve a blue- 
and redshifted peak only at 3.4mm while at 1.4mm we only 
resolve a blueshifted peak. The peaks at both wavelengths are 
resolved, and the source MM2 is the most Ukely driving source. 
Its central spectrum (Fig.|4]3) shows clear blue- and a redshifted 
peaks. This outflow also seems to be aligned close to the l.o.s., 
according to the position of the blue- and redshifted emission 
peaks without any emission at the Vlsr- 

The outflow SIO-c is detected in both SiO transitions, al- 
though at 3.4mm the redshifted lobe is faint. The outflow is 
aligned in the east-west direction, its axis defined as the straight 
line joining the blue and red emission peaks at 3.4mm. The 
source MM4 lies on the outflow axis and between the lobes, thus 
this may be the driving source. The SiO (5^) spectrum from 
this outflow (Fig. |4j;) shows broad blue- and redshifted peaks. 
The blue peak is broad (~ 5.5 km s"') and asymmetric towards 
bluer velocities, its maximum at -44 km s" ' , while the redshifted 
peak is symmetric and narrower, with a FWHM linewidth of 
4.8 km s"' peaking at -36kms"'. 

For SlO-d we do not detect any millimetric source that 
could be driving it. This outflow is best traced at 3.4mm in 
the SiO (2-1) transition, where both blue- and redshifted lobes 
are clearly seen, while at 1.4 mm just the redshifted emission is 
traced. The central spectrum for this outflow (Fig. |4}l) shows 
broad (~ 5.5 km s"') double-pea ked emission centered at the 
systemic velocity. On larger scales lRidge & Moorel(l200H) found 
a CO outflow spanning 1 .73 pc in the northwest-southeast direc- 
tion, the same as SlO-d (see Section|4j2]i. 

SlO-e is even more intriguing because it is the biggest of the 
five outflows and we cannot identify a driving source for it in 
the millimetric or the IR. Like SlO-d it lies in the northwest- 
southeast direction, matching the large-scale CO emission. Its 
SiO (2-1) spectrum shows broad blueshifted emission, peaking 
at -40.4 km s ' with a FWHM of 5.6 kms '. The redshifted 
peak is narrower with a FWHM of 3.4km s"' and with its max- 
imum at -38 kms"'. Its driving source should be located rela- 
tively far from the two main sources MMl and MM2, at ~10'* 
AU from them. 



3.4. Sulfur Dioxide emission 

Figure |5] panels (a) and (c) show the integrated Sulfur Dioxide 
(SO2) emission detected towards W3IRS5 at 1.4 mm and 
3.4 mm respectively, while in panels (b) and (d) are the velocity 
maps for each transition. The SO2 emission envelops the whole 
region where MMl, MM2, MM3 and MM4 are detected, which 
can be clearly seen in Figure|5j;. We see a clear velocity gradient 
in the southeast-northwest direction, indicating that the enve- 
lope of W3 IRS5 rotates as a whole and is likely gravitationally 
bound. This is more noticeable in the SO2 (83,5-92,8) velocity 
map (Fig. |5]i), where the velocity gradient goes smoothly from 
-42 to -37 kms"'. 

The peaks of integrated SO2 emission at 1.4 mm (Fig. |5^) 
He close to the sources MMl, MM2 a nd MM4. As SO2 is also 
an indicator of shocked material (e.g. IHelmich et al.lll994b this 
strengthens the argument for the detection of the outflows traced 
by SiO, and supports the statement in Section 13.31 that those 
sources drive the outflows SlO-a, SlO-b and SIO-c, respectively. 

The SO2 (222,20-22 1,21) velocity map (Fig. [S}?) shows a 
strong jump of ~ 4 kms"' in a thin spatial region southeast 
of the two main millimetric peaks. The velocity shift is from 
35 kms"' to ~ -39 kms"' -the systemic velocity- in a nar- 
row strip of ~ 1.1" X 0.2". The region north of the jump is 
where most of the gas is blueshifted while the region south of the 
jump is completely redshifted. This "velocity jump" (see Sec. 
14.31 ) may be the signature of two molecular flows moving in op- 
posite direction, rammin g each other and compressing the gas in 
the region of the ju mp (Kles sen et al. I I2OO5I: iHeitschet al.ll20()l 
l2006l:lPeretto et al.il2006. ,20071) ! 

4. Discussion 

4.1. Continuum Sources 

With L = 2 X lO^Lo, W3 IRS5 is a high-mass star-forming re- 
gion, but the (proto)stellar content in its core is unknown be- 
cause of high obscuration. Since it is a hypercompact Hii region, 
it already has a large (proto)stellar component. While it is not 
possible to directly measure the masses of the (proto)stars with 
millimeter data, it is possible to derive the masses of the cir- 
cumstellar structure. The low masses that we calculate from the 
1.4 mm data can be attributed to a circumstellar structure, while 
the surrounding envelope likely contributes to the larger masses 
derived from the 3.4 mm data. However, we must take into ac- 
count the "short-spacing problem" when calculating the masses 
of the continuum sources. This caveat, inherent to interferome- 
ters, means that a high percentage of the flux is filtered out and 
lost. This eff'ect is more severe in the extended configurations 
(long baselines), as in our case. 

Single-dish submil limetric measurem ents for W3 IRS5 at 
850//m with SCUBA ( Moore et alj|200 7l) indicate a peak inten- 
sity Iv ~ 9.6 Jy beam"' with a beam of 14". Considering a de- 
pendence ly oc for the peak intensity, that value corresponds 
then to \y ~ 1.4 Jy beam"' at 1.4 mm and ly ~ 36mJybeam"' 
at 3.4mm. In our 1.4mm continuum map the flux density re- 
covered within the SCUBA beam is Sv ~ 0.12Jy taking the 
average contribution of the free-free emission into account (see 
Table |3]l, meaning that ~ 91% of the flux is lost. Therefore the 
calculated masses at 1.4 mm are actually a lower limit for the 
current mass of each source. On the other hand, at 3.4 mm we 
recover S^ ~ 36 mJy, so with our calibration accuracy, we do not 
miss too much of the flux within the SCUBA beam. However, 
within the PdBI main beam of ~ 55" at 3.4 mm we only re- 
cover ~ 73 mJy, while iMoore et al.l (l2007h reports an integrated 
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flux Sy ~ 200 Jy, implying that although we are recovering most 
of the central emission, we are filtering out the more extended 
emission. Nevertheless, at 3.4 mm the obtained values for the 
masses are somewhat more representative of the mass of their 
envelopes. 

The SO2 velocity signature across the core is indicative 
of overall rotation of the gas associated with MMl to MM3. 
Assuming equilibrium between the gravitational and rotational 
force at the outer radius of the cloud: 



Table 4. Millimetric, radio, NIR and MIR counterparts. 



5v r 



M„,[M0] = 1.13 10"^ X 5v2[kms"'] X r[AU] 



(1) 



(2) 



where r is the cloud radius (~ 2" - 4000 AU) and 6v is half the 
velocity regime observed in the first-moment map (~ 2.5 km s"' ; 
Fig. El). 

The gravitationally supported mass implied by this observed 
rotation is Mrot ~ 30 M©, which is on the same order as the de- 
rived masses from the 3.4 mm continuum emission. Hence, these 
sources are likely gravitationally bound and may still be accret- 
ing from an envelope not seen in our data, further increasing 
their masses. Also, from the distribution of the gas it is likely 
that sources MMl, MM2 and MM3 are sharing a common en- 
velope, of which we are tracing the inner parts at 3.4mm. This 
envelope would have much more of the mass than in the imme- 
diate region surrounding each source, but it is being filtered out 
by the interferometer 

The different percentage of flux filtered out in each wave- 
length band can be explained the same way as the big (over a 
factor ~ 10) difference between the calculated masses at 1.4 mm 
and 3.4mm. Both wavelengths were observed at the same time 
with the same interferometer configurations, meaning that al- 
though the ground baselines are the same, the uv coverage mea- 
sured in units of wavelength (k/i) at 3.4mm is more compact 
than at 1 .4 mm, therefore tracing more extended components in 
the region. 

Even if we are recovering most of the compact flux at 
3.4 mm, we are likely mapping the surrounding envelopes close 
to the already formed (proto)stellar objects and therefore can- 
not state the actual masses of the (proto)stellar objects but only 
give a lower limit of the mass of the gas and dust associated 
largely with the individual sources. The same can be applied for 
the column density, however in this case -despite the caveats 
mentioned- the values we obtain are more realistic for a region 
like W3IRS5. Our results of N{H 2) ~ lO^^cm^^ at 3.4mm, 
are 2-3 magnitudes higher than for low-mass star-forming re- 
gions, for which the derived column de nsities at ~3 mm are of 
flie or der N(H2) ~ 10^2 -23 ^ . [Hariunpaeae & Mattilal 

Il996t iMotte et al.lll998t iKontinen et aimodi . Also our results 
are about two orders of magnitude higher than the critical col- 
umn density of 1 g cm"^ required to for m high-mass stars un- 
der M ilky Way conditions calculated by iKrumholz & McKeg 
( I2OO 8). which corresponds to N{H2) - 3 x 10^^ cm"^, suggest- 
ing that within the sources we are mapping, (proto)stars of over 
IOM0 are forming. 

Radio observations for this region are also available. 
Comparing the positions of the five radio sources detected 
with the VLA at 43 GHz (g-band) bv Ivan der Tak et alj (l2005h 
and the positions in Table |2] we can identify the following: 
MMl ^ Q5, MM2 ^ Q3 and MM6 ^ Q4 (see Fig. EJ. The 
source MM3 lies between Ql and Q2 and is not a compact 



MM 



Ql 



NIR" MIR" 



5 

3 

1+2 

1+2 + 3 
4 



1 
2 

2a 



" From'van der Tak et al.' (200^ 
* From .Megeath et al., L2005.) 



source but rather stretched out, so it appears to be the joint con- 
tribution from Ql and Q2, which could not be individually re- 
solved with our spatial res olution (see Table[4|. 

At shorter wavelengths. lvan der Tak et al. I (l2005l) found three 
MIR sources in W3 IRS5 with the Keck I telescope, while 
Megeath et al. (2005) found seven NIR sources with the Hubble 
Space Telescope, three of them which match the aforementioned 
MIR sources. 

Because absolute astrometry was not available for these NIR 
and MIR sources, to compare them with our MM sources we 
refer to their relative positions, shown in Table|5] The separation 
between MMl and MM2 is 1.16" + 0.10" at a position angle 
(RA.) of 214 + 5 degrees, being the only match for the pair MIRl 
and MIR2 with a mean separation of 1.12" ± 0.07" and a RA. 
of 216.8 + 1.7 degrees, and for the pair NIR 1 and NIR 2 with 
a separation of 1.23" and a RA. of 217 degrees. Similarly, the 
pair MM1-MM6 with a separation of 2.65" + 0.10" and F.A. = 
187 + 2 degrees is the only match for the pair MIR1-MIR3 with 
a mean separation of 2.70" ± 0.05" and a RA. of 187.2 ± 0.7 
degrees and for the pair NIR 1-NIR 3 with a separation of 2.7" 
and P.A. = 189 degrees. Therefore we can identify MIRl and 
NIR 1 with MMl, MIR2 and NIR 2 with MM2, and MIR3 and 
NIR 3 with MM6. Comparing the separation of 0.94" ± 0.10" 
and P.A. = 235 between MMl and MM3, with the separation 
of 0.98" and P.A. = 224 between NIR 1 and NIR 2a, we can 
also identify source NIR 2a with MM3. For the rest of the NIR 
sources we do not find a millimetric counterpart (see Table|4]l. 

The only millimetric source without any NIR, MIR or radio 
counterpart is MM4, hence it is a new detection. Its mass could 
not be corrected by the free-free emission contribution and it is 
approximately half as massive as MM3. 

Considering the region encompassed by MMl and MM2, 
we estimate a lower limit for the (proto)stellar density. 
Assuming spherical symmetry, we detect the four sources 
MMl, MM2, Ql and Q2 (joined in MM3) within a 1.16" or 
~2100AU diameter region, corresponding to a (proto)stellar 
density of ~ 7 x lO^protostarspc"^. Even counting just the 
three resolved sources MMl, MM2 and MM3 we obtain 
~ 5 X 10^ protostarspc"^. It is not the first time that (proto)stellar 
densities higher than the typical 10"* stars pc~^ for you ng clusters 
(iLada & Ladall2003h are found. ' Beuflier et all (l2007h calculated 
a density of ~ 1.4 x 10^ protostars pc~^ in the system G29.96- 
0.02. Although high, this (proto)stellar density is still below the 
"critical" val ue of ~ 10*^ protostar s pc"^ predicted by the merg- 
ing scenario (iBonnell et al.lll998l) . However it is the first time 
that densities higher than the required ~ 10''protostarsp c~^ to 
induce high-mass close binary mergers (Bonnell & Bate l 2005|) 
are obtained. While the observations of W3 IRS5 show no evi- 
dence of the presence of such binary systems, we have shown 
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that such densities can be achieved in high-mass star-forming 
regions, and with further improvement in the spatial resolution, 
(proto)stellar densities as high -or even higher- may be a com- 
mon phenomenon. 



Table 5. Relative astrometry of the known sources in W3 IRS5. 





Sep 


PA 


Sources 


(arcsec) 


(deg) 


MM1-MM2 


1.16 


214 


MIR1-MIR2;;. 


1.12 


216.8 


NIR1-NIR2^ 


1.23 


217 


MM1-MM6 


2.65 


187 


MIR1-MIR3 


2.70 


187.2 


NIR1-NIR3 


2.7 


189 


MM1-MM3 


0.94 


235 


NIRl-NIR2a 


0.98 


224 



" Mid-infrared sources from lvan der Tak et all j2005h 
* Near-infrared sources form lMegeath et al] ( l2005h 



W3 IRS5 is surrounded b y a dense embedded cluster of 
intermediate- to low-mass stars: iMegeathetaildlQgd ) calculated 
a stellar surface density of ~ 0.25 x 10'*pc~^ for a cluster ra- 
dius of 21", a quarte r of the ~ 10"* pc"^ density calculated by 
iMegeath et alJ (l2005h from the NIR sources they found with a 
maximum projected separation of ~ 5,600AU. Again within a 
~ 2100 AU diameter region, we calculate a (proto)stellar surface 
density of ~ 4 x lO'^protostarspc"^ considering just the three 
millimetric sources previously mentioned, revealing a density 
gradient towards the center of the cluster Assuming that this 
embedded cluster around W3 IRS5 forms an independent sub- 
cluster within the W3-Main region, the fact that there are already 
low-mass stars in the outskirts of W3 IRS5 while high-mass stars 
are still forming in its center supports the hypothesis that the 
low-mass stars for m before their high-mass counterparts (e.g., 
iKumar et alj|2006l) . 

The projected separation between sources ranges from 
~750AU for MM2-MM3 to ~4700AU for MM1-MM6. 
Compared to the median radius of 40,000 AU for the farthest 
outlyin g member of the 14 trapezia identified bv lAbt & Corballvl 
(I2OOOI) . W3 IRS5 has significantly smaller projected separa- 
tions than those observed in optically visible trapezia. Therefore 
we have several compact sources sharing a common envelope 
from which gas may still be accreting in a region with a high 
(proto)stellar density beyond the point where close binary merg- 
ers can be induced, and with column densities higher than the 
theorized threshold to form high-mass stars. Thus, conditions are 
favorable for the formation of high-mass stars. 



4.2. Outflows 

For the five SiO outflows mapped in this region, we identify the 
driving sources of just three of them. Outflows SlO-a, SlO-b and 
SIO-c are driven by sources MMl, MM2 and MM4 respectively, 
while for outflows SlO-d and SlO-e we do not d etect candidates 
for th eir driving sources. For SiO-e however, 'Tieftru nk et al.l 
([1998) have detected a NH3 (1,1) emission peak (labeled 1 in 
their Figures 4 and 5) near the center of this outflow, which may 
be hosting the driving source of SiO-e. 



The determination of the outflow parameters is subject to a 
big number of errors, including that it is often difficult to sepa- 
rate the outflowing gas from the ambient gas and that the incli- 
nation angles of the flows are also often unknown. In this case 
we also have the caveat that high amounts of flux are filtered out 
by the interferometer. To get an estimation of the ef fect of this 
we compare our results with the measurements of iGibb et al.l 
(2007). They report an integrated intensity of 2.4Kkms^' for 
SiO (5^), while we recover only 0.4Kkms ' within our pri- 
ma ry beam. A l thoug h apparently the integrated intensity given 
by iGibb et aTl (|2007|) covers a slightly larger region, the main 
emission stems from an area covered by approximately twice 
our 22" primary beam at 217 GHz, implying a great fraction of 
missing flux also in SiO (5^). On the other hand, for SiO (2-1) 
they measure 2.5 Kkm s"' but with most of the emission within 
a region covered by our 55" primary beam at 87 GHz, in which 
we measure 1 .7 K km s" ' . Again, we are losing a high amount of 
flux, although less than in SiO (5-4). 

We derived the physical parameters for the small-scale out- 
flows detected here. To obtain the H2 column densities we first 
have to calculate the SiO column densities for each outflow, for 
which we follow the LTE approach described by llrvine et akl 
([T987h . Since we do not have estimates of the excitation tem- 
perature, we ass ume a usual value for the outflows of 30 K for 
our calculations jBeuther et al.ll2002bl) . resulting in the formula: 



A^s,o(2-i) = 1.3x10'^ rdv 



(3) 



where J Tdv is the SiO (2-1) integrated intensity in Kkms"'. 
We used the SiO (2-1) data because it has recovered the most 
flux. The obtained results are shown in Table |6] Adopting a 
SiO abundance ratio to obtain the H2 column density is not 
straightforward. The SiO abundance has a big un certainty, with 
repor t ed values ranging f rom SiO/H2 ~ 10"^ dMikami et al.l 



1 19921 IZhangetal.1 Il995h for shocked energetic regions to 
SiO/H, ~ 10-'^ (IZiurvs&Friberg|ll987t IZiurvs et al.lfl989l) in 
quiescent, cold dark clouds. To take into account such uncer- 
tainty, our calculations were done assuming SiO/H2 = 10"^ 
and SiO/H2 = 10"^, a range suitable for a high-luminosity 
region such as W3IRS5. Although the values we obtain are 
too low for a region Uke W3 IRS5 with a total luminosity of 
~ lO^Lo and are more typical in regions with lower luminos- 
ity (see e.g. Hidge & Moore 2001; Gibb et al. 2007; Palau et al] 
I2OO7I; iBeltran et al.ll2008h . we have to keep in mind that aside 
from the caveats mentioned above, we are tracing smaller scales 
opposed to large-scale outflow. 

In general, our parameters are several orders of magnitude 
lower than for the CO outflow detected by Ridge & Mo or3 
(I2OOII) even assuming extremely low SiO abundances. From 
the CO measurements, the outflow has a mechanical en- 
ergy Eco ~ 10"*^ erg, while for a SiO abundance of 10"^, 
adding the contribution of all the outflows we obtain that 
Esio/Eco ~ 4 X 10 From the values in Table |6] it is evident 
that the rest of the parameters follow the same trend. Only at 
unrealistically low SiO abundances are our results comparable 

with the CO outflow. 

The SiO outflow that lGibbet all (12007) detect is aligned in 
the east-west direction, and the CO outflow also has an east- 
west component peaking at ~35" east of the main source ex- 
tending to the west, as well as a NW-SE component. In Fig. |3] 
we can see that the SiO (2-1) emission follows a similar distri- 
bution, with east-west (outflow SIO-c) and NW-SE (outflows 
SlO-d and -e) components. As mentioned before we do not de- 
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Table 6. Properties of the small-scale outflow features 



Outflow M<„„ E M„„, 

(Me) (erg) (Mgyr'') (Mg yr'' km s'') (L^) 

SiO/H2 = 10"^ 

SlO-a + SlO-b^ 004 4(43) 4(^5) 2(5) 0^ 

SIO-c 0.03 5(43) 5(-5) 7(-4) 0.80 

SlO-d 0.03 1(43) 3(-5) 2(-4) 0.09 

SlO-e 0.08 3(43) l(-5) 6(-5) 0.03 
SiO/Hj = lO-** 

SlO-a + SIO-b° 4 4(45) 4(3) 2(3) 25~ 

SIO-c 3 5(45) 5(-3) 7(-2) 80 

SlO-d 3 1(45) 3(-3) 2(-2) 9 

SlO-e 8 3(45) l(-3) 6(-3) 3 



" Treated jointly because they cannot be differentiated in the red- 
shifted SiO (2-1) emission. Also, because they are aligned close to the 
1.0. s., the parameters for which the outflow size is needed were calcu- 
lated assuming an average size based on the other three outflows. 

tect driving sources for SlO-d and -e. Because of their orien- 
tation these two outflows could be delineating the high-density 
regions of the large-scale CO outflow and although we cannot 
exclude this option, we believe that they are two separate fea- 
tures. Furthermore, Tieftrunk et al. (1998) detects an ammonia 
signature close to the projected center of SlO-e which may be 
hosting its driving source. With our data we cannot confidently 
point to any of our detected dust sources as the driving source 
of the large-scale CO outflow and it is possible that this CO out- 
flow is the result of the mix of all the small-scale outflows and 
therefore might have driving sources opposed to a single driving 
source. 

As shown in Section [33] the outflows SlO-a and SlO-b are 
aligned very close to the l.o.s.. This situation explains the fact 
that despite the high column densities toward MMl and MM2 
(see Table |2]i, which correspond to visual extinctions of several 
thousand magnitudes, NIR and MIR counterparts are detected. 
The cavities carved out by the outflows have reduced the extinc- 
tion along the l.o.s., allowing us to see hot dust in the interior of 
the protostellar envelopes in the NIR and MIR, either directly or 
through scattered light. 

In the case of SiO-b, the outflow also would explain the oc- 
currence and location of NIR 2b, below NIR 2/MM2. That near- 
infrared source can be interpreted as scattered light from the out- 
flow cavities. Similarly in the case of source NIR 2 a, it is close to 
the pr oper motion radio source IRS 5a detected bv I Wilson et alj 
(12003) with the VLA. They propose that IRS 5a/NIR 2a was 
formed by or in an outflow originating south of its position, 
which agrees with the localization of our outflow SiO-b. 

From the H2O masers they detect, llmai et aL I dlOOQ) suggest 
two outflows in the north-south direction in this region. However 
from our data we cannot identify counterparts for the position of 
the H2O masers in this region, nor outflows in a north-south 
direction. 

4.3. Velocity jump 

Figure |6] shows a postion-velocity diagram across the velocity 
jump seen in SO2 (222,20-22 1,21). Two distinct velocity compo- 
nents are seen, one peaking at 35 km s ' and the other peak- 
ing at 39kms"^ spatially located at both sides of the veloc- 
ity jump, denoted by the black dashed line. The velocity jump 
has also been confirmed at higher frequencies, in SO2 (8-7) and 



HCN (4-3) observations taken with the SMA towards W3 IRS5 
(T. Bourke, priv. comm.), and could be explained by the the- 
ory of gravoturbulent star formation (H eitsch et al. 2005, 200^ 
IKlessen et al J 120051: iMac Low & Klessenli2004i and references 
therein), which states that protostellar cores might form by tur- 
bulent ram pressure compression in the regions behind the fronts 
of converging turbulent flows. 

From Figures |2] and I5J5 we can see that the main sources of 
W3 IRS5 are located next to the region where the velocity jump 
occurs, distributed roughly parallel to the jump. This is what 
would be expected if their collapse were triggered by the com- 
pression of co nverging flows. Furthermore, numerical models 
(Kl essen et al. 2005) predict that one of the observational signa- 
tures of such compression by ram pressure is localized maxima 
of the l.o.s. velocity dispersion in the low column density gas in 
the outskirts of the core. According to the continuum maps (Fig. 
|2]l, the region where this velocity jump occurs is in the outskirts 
of the cores, where the dust emission and the gas column density 
are low. This suggests that we may be looking at two molecular 
flows converging and triggering the star formation in W3 IRS5. 

At larg er spati al scal es, a similar signature was found by 
iPeretto et al.1 (120061 l2007h in their study of NGC 2264-C at mil- 
limetric wavelengths, where they found in N2H^ and H'^CO^ 
emission a velocity jump of ~ 2kms"' over a very narrow 
spatial region. They modeled this jump as the signature of the 
large-scale, axial collapse of NGC 2264-C along its long axis 
and towards its center, with an inclination angle of ~ 45°. It 
would also trace a possible dynamical interaction between pro- 
tostellar sources at its center, where they detected a millime- 
ter continuum peak. The case of W3 IRS5 could be similar 
but observed with higher spatial resolution. Making an anal- 
ogy between iPeretto et al.l (2007) and W3 IRS5, sources MMl, 
MM2 and MM3 would correspond to their sources CMM3 and 
CMM13, and we would just be mapping the region of the jump 
from redshifted emission to the systemic velocity, which corre- 
sponds in their case to the transition from CMM4 to CMM13- 
CMM3. 

5. Summary and Conclusions 

In the continuum maps of W3 IRS5 obtained with the PdBI at 
1.4 mm and 3.4 mm we detect 5 individual sources. Their calcu- 
lated absolute masses are strongly affected by the spatial filter- 
ing inherent to the interferometric technique, however from the 
relative separations of the three strongest sources we propose 
the scenario of a Trapezium-like system, with at least 3 objects 
enclosed within a ~2000AU radius, with a high (proto)stellar 
density of ~ 5 x lO^protostarspc^^*, and sharing a common en- 
velope from which they may still be accreting gas. 

Four of the five millimetric sources are identified as coun- 
terparts of sources previou sly detected in N IR, MIR and ra- 
dio wa velengths by Megeat h"et al.l (l2005h and Ivan der Tak et"al] 
(I2OO5I) . The remaining fifth source (labeled MM4) is a new de- 
tection. 

Thanks to the spatial resolution of ~0.36" we can for the 
first time trace the small-scale flow, detecting five molecular out- 
flows, identifying the driving source of three of them. Although 
we find on short spatial scales the same alignment and axes as the 
outflows on large spatial scales previously detected in SiO and 
CO, within our accuracy and the big uncertainties introduced, 
we do not find strong evidence to support that we are tracing the 
inner and denser regions of those outflows, but we cannot rule 
out that possibility. However these small-scale outflows appear 
to be different features. 
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Two of the SiO outflows are close to the l.o.s. direction and 
each one is being driven by each of the stronger millimetric 
sources detected. This configuration would explain why these 
two objects are detected from NIR to radio wavelengths because 
cavities being carved out by outflows nearly along the l.o.s. have 
lower extinction, allowing the detection in the NIR and MIR of 
the hot dust near the collapsing protostar, while the long wave- 
length emission traces the dusty envelope that still surrounds it. 

The SO2 emission and velocity structure indicate that 
W3 IRS5 still has a gravitationally bound envelope. Also in SO2 
we detect a velocity jump of several kms ' in a very narrow 
region adjacent to the continuum sources. This could be the sig- 
nature of two converging molecular flows, compressing the gas 
and triggering the star formation, as predicted by the theory of 
gravoturbulent star formation. 
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RA (J2000) 

Fig. 1. Continuu m image (greyscale) of the W3 region at 4.9 
GHz adapted from lTieftrunk et al] (1 19971) . Contours trace the ob- 
served PdBI continuum emission at 3.4 mm, starting at the 3cr 
level in Icr steps. Shown are the Hii regions W3 B (compact), 
W3 F (ultracompact) and W3 M/W3 IRS5 (hypercompact). W3 
B lies at the border of our primary beam, thus the contours above 
and below W3 B are artifacts of the deconvolution process, and 
do not have any physical meaning. 
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Fig. 2. Left and center panels: Continuum maps of W3 IRS5 at 
1.4mm and 3.4mm respectively. Contour levels for the 1.4mm 
continuum start at 4cr increasing by 2cr steps, and at 3cr increas- 
ing by Icr steps for the 3.4mm continuum. Beams are shown 
at the bottom left of each panel. Filled triangles mark the milli- 
metric sources detected (see Table |2]i. Right panel: Relative po- 
sitions of the different sources in W3IRS5. Filled triangles are 
the milli metric sources detecte d, stars are the NIR sources de- 
tected by Megeath et al.'(2005') and crosses are the radio sources 
detected by van der Tak et al. (2005). The grey square marks 
source MM5, the joint contribution of sources MM2 and MM3. 
In Table ID are the identified counterparts, discussed in the text. 



Fig. 3. Silicon Monoxide emission detected towards W3IRS5. 
The left panel shows the 3.4mm SiO (2-1) transition and the 
right panel shows the 1 .4 mm SiO (5^) transition. Solid con- 
tours outline the blueshifted emission and dashed contours the 
redshifted emission. The blueshifted emission ranges from -58 
to -43kms"' and the redshifted emission ranges from -36 to 
-28 kms"'. In the left panel contouring starts at the 4cr level 
and in the right panel at the 3cr level, increasing in Icr steps 
for both. The beams are represented at the bottom right corner 
of each panel, respectively. Labeled SlO-a through -e are the 
five outflows detected, filled triangles are the millimetric sources 
shown in this paper and the stars mark the NIR sources from 
iMegeafli et all ( l2005l) . 
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Fig. 4. SiO spectra of the five detected outflows. In panels (a), 
(b) and (c) are the SiO (5^) spectra of outflows SlO-a, -b and -c 
respectively. In panels (d) and (e) are the SiO (2-1) spectra of the 
outflows SlO-d and e respectively and in panel (f) is the spectrum 
of the SiO (5-4) integrated emission in a 4" x4" region centered 
at the phase center The vertical dotted line in each panel marks 
the systemic velocity of -39kms"'. In panels (c) and (e) the 
solid line is the blueshifted emission and the dot-dashed line is 
the redshifted emission. In panel (f) the long-dashed lines de- 
mark the velocity regime for the ambient gas. The spectra were 
taken toward the center of the outflow for SlO-a, -b and -d, and 
toward the wings for SIO-c and -e. 
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Fig. 5. Sulfur Dioxide emission detected towards W3IRS5. At 
1.4mm is the SO2(222,20-22i.2i) transition (upper half) and at 
3.4 mm is the 802(83 5-92,8) transition (lower half). Panels (a) 
and (c) show the integrated emission, with contour steps of 
0.2 Jy beam"' km s"' for (a) and O.lJybeam"' kms"' for (c). In 
panels (b) and (d) is the velocity distribution of the emission, 
from -48 to -33kms"' for (b) and from -42 to -37kms"' for 
(d), in both panels the contour step is 0.5 km s ' . The respective 
beams are shown at the bottom-right corner of each panel, re- 
spectively. Filled black and white triangles are the mm sources 
from Table |2l The dashed white line in panel (b) marks the ve- 
locity jump discussed in the text. 



T 




offset [arcsec] 

Fig. 6. SO2 (222,20-22 121) position velocity diagram along the 
line shown in Fig. ^h) with a white arrow. The velocity jump 
described in the text can be seen here, where an emission feature 
is peaking at ~ -35kms"' while another feature is peaking at 
~ -39kms"', spatially at either side of the jump, denoted by 
the black dashed line. Contours start at 20% of the maximum 
flux and continue in 10% steps. 



